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Q ■ Abstract 

We performed an analysis on the detection of a long-lived slepton at a linear 
collider with yfs = 500 GeV. In GMSB models a long-lived NLSP is predicted 
for large value of the supersymmetry breaking scale \fF. Furthermore in a 
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lived staus with masses up to the kinematical limit of the machine. We also 
present our estimation of the sensitivity to the stau lifetime. 
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large portion of the parameter space this particle is a stau. Such heavy charged 
particles will leave a track in the tracking volume and hit the muonic detector. 
In order to disentangle this signal from the muon background, we explore 
kinematics and particle identification tools: time of flight device, dE/dX and 
Cerenkov devices. We show that a linear collider will be able to detect long- 
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I. INTRODUCTION 



In many supersymmetric models charged long-lived particle can exist. In models with 
R— parity violation, this happens when the lightest supersymmetric particle (LSP) is a 
slepton and the R— parity violation term is small |l| . Another class consist of models where 
the LSP is not the usual (like in mSUGRA) U(1)y gaugino but rather a SU{2) gaugino; in 
this case the next-to-lightest supersymmetric particle (NLSP) is the lightest chargino, almost 
degenerate in mass with the neutralino LSP, and will decay after traveling centimeters or 
meters into the LSP plus a soft lepton or pion ||. This type of model arises rather naturally 
in the anomaly mediated SUSY breaking scenario ||, however in this scenario the mass 
splitting is typically larger then the pion mass, leading to decay length of the order of 
ct ~ 10 cm at most. Another class of models where long-lived charged particles occur 
is in gauge mediated supersymmetry breaking (GMSB) models |4]||. Depending on the 
value assumed for the SUSY breaking scale, V~F, the gravitino is the LSP. Moreover, due 
to the weak coupling of the gravitino to the SUSY fields, the NLSP could in principle be a 
long-lived massive particle. Throughout this paper we will concentrate on the detection of 
long-lived charged particles in the GMSB framework. The techniques presented, however, 
are applicable to general non-strongly interacting heavy long-lived charged particles []. 

In GMSB an intermediate sector is responsible for communicating Supersymmetry break- 
ing to the MSSM sector. In these models the gravitino mass is given by, 
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M r = ^ = * J 2.37 eV (1) 



where M p i is the reduced Planck mass. For the yF ~ 100 TeV (to be compaired with 
10 11 GeV, a typical mSUGRA value of SUSY breaking scale), Mq ~ eV, making the phe- 
nomenology quite interesting. 

Moreover, because of the very weak interaction of the gravitino, all supersymmetric 
particles will decay into the NLSP and this will decay into the gravitino (which will escape 
detection) and its standard model partner. Thus, the nature of the NLSP and its decay 
length will play a fundamental role in the phenomenology of the model. Because soft terms in 
GMSB models are generated by the gauge couplings, the NLSP is usually either a neutralino 
or a stau. Many phenomenological studies have been made for Tevatron and LEP in the 
context of this model [0-|9j . More recently some studies for LHC and Run II at Tevatron also 



have been addressed fioHT2f . For the next linear collider (NLC) the case for a neutralino 



NLSP, both long- and short-lived, was considered by Ambrosanio [13| and the case for a 
stau NLSP with prompt decay was considered by Kanaya |L4j]. 

In this work we improve over our previous analyses P^JI3[] in the search for a long-lived 
stau at NLC, presenting here our estimations for a lifetime measurement. In the GMSB the 
NLSP lifetime is given by, 



1 In the Ref. ||, the authors raise the issue of hadronization effects for color-triplet particles. Some 
of them will hadronize into neutral exotic mesons. Besides that, the inelastic hadronic reactions 
can change the charge of the mesons. 
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or = 16,- [-f—l (H^YlOkm. (2) 
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From this expression, one can see that the measurement of the NLSP mass and lifetime 
will determine the value of the fundamental SUSY breaking scale v~F. From perturbative 
arguments it is possible to set a lower limit in \J~F for the complete set of parameters in 
the GMSB framework. Nevertheless, for the range of stau masses under consideration, it 
is always possible to find parameters which gives a cr ~ fim. On the other hand, there is 
no solid theoretical argument for an upper limit on y/~F. However, a LSP gravitino with 
mass higher than few KeV is disfavored in some cosmological scenarios for over-closing the 
Universe PI . This translates to roughly ^/F < 3000 TeV. 



The rest of this paper is organized as follows. In section II we describe our simulation 



using the program ISAJET |L7j to generate stau pairs with the effects of a time of flight 
device, dE/dX and a Cerenkov device to detect a heavy particle. In section III we use the 
stau pair production process to extract limits on the lifetime of stau which leads to limits 
on the SUSY breaking scale in a GMSB model. In section IV we present our conclusions. 



II. SELECTION CRITERIA FOR STAU PAIR PRODUCTION 

Stau pair production at a linear collider provides a clean search environment. Further- 
more, the production cross section is model-independent, depending only on the mass of 
the staus and on the mixing angle between the left- and right-handed superpartners. In 
Fig. |T| we show the pair production cross section (normalized to <x w = 450 fb) for the left- 
and right-handed states as a function of stau mass. We can see that the stau pair cross 
section is smaller than due to its scalar nature and it rapidly drops when we approach 
the kinematical limits of the accelerator. Nevertheless, we note that for masses around 240 
GeV we still have cross section of C(10) fb, which should be observable provided that the 
background is manageable. 

We will first describe the selection criteria without particle identification which will be 
used in later sections to extend the range of sensitivity to the full beam energy. The signal 
we are looking for is a back-to-back tracks with corresponding hits in the muon chamber. 
With this requirement tracks from 7t,K,p and e are removed. To reject the two photon 
process of 77 — > we note that the pair in this process tends to have a low 

invariant mass and be boosted along the beam pipe. Thus, we require the following cuts: 

1. cos 9 < 0.8, to guarantee good track quality; 

2. \p\ > 0.5E^ eam and 

3. \pt ot \ < 0.25£ beam . 

After these cuts, the two photon initiated muon pair production is estimated to be 1.4 fb. 
We are then left with muon pair production e + e~ — > as the main source of background. 

In order to reduce the muon pair background we shall explore the heavy mass of the 
staus. In a e + e~ collider the energy in the center-of-mass is fixed, so in a pair production 
process the energy of the final particles is also known. It is well known, however, that 
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in a high energy linear collider beamsstrahlung and initial state radiation effects become 
important and the effective energy of the reaction is not fixed but presents a spectrum. 
In ISAJET, these effects have been implemented using the parameterization described in 



Ref. [jig] , where the radiation spectrum is well approximated by one photon emission from 
one of the initial state e . Assuming also that the direction of the emission is along the 
beam line, the mass estimate for each track in the laboratory frame is given by: 



M 2 = ^- + (3p z - \p\ 2 , (3) 




S = a (l-|A|), (4) 
/3 = A/(2- |A|) , (5) 

where A = p\ /-E'beam is the net momentum in the beam line direction, (3 is the boost 
parameter and \/I is the center-of-mass energy of the two tracks. 

In Fig. is shown the cross section distribution as a function of M 2 which is estimated 
according to Eq. (|3]). In this plot we can see the muon distribution peaking at zero mass 
with a tail from beamsstrahlung. We also see the distribution of staus production for several 
values of masses. The momentum resolution is taken to be Spt/pt = 5 x lO~ 5 p T (GeV). 
Based in this plot we use the following cut, 

4. \M 2 -M?\< 3000 GeV 2 , 

where M f is a variable parameter in the search. The resulting efficiency after this cut is 
shown as the dotted line of Fig. ^|a. This is our basic strategy to reduce the muon background. 
To further improve the sensitivity, we study particle identification. 



A. Time of Flight 

A time of flight (TOF) device can be used to identify heavy tracks. In our study we 
considered a linear collider with 1.4 ns of bunch separation. In the large detector scenario 
(r = 2 m) the mean time of flight for a massless ((3 = 1) particle is around 6.7 ns. Assuming 
that we do not know which bunch crossing a given event is coming from, every 1.4 ns of 
TOF delay is consistent with the massless muon. We simulated the effect of a 50 ps error 
in the time of flight measurement, applying the following cut, 

5. At > 0.13 ns, 

where At is the time of flight difference between a (3 = 1 and a massive particle, modulo 1.4 
ns. 

This cut corresponds to about 2.5a, so that about 1% of the muon background are kept. 
Applying this cut to both tracks we can relax cuts 2 and 3 which extends the mass range 
to the full beam energy. The efficiency of this cut as a function of mass is shown in Fig. ^a 
as the dashed line. The corresponding sensitivity is shown in Fig. 0b. We see that for some 
values of the mass we have lower efficiency when the time delay with respect to (3 — 1 case 
becomes equal to the bunch spacing. However, the efficiency is never below 0.3 because the 
actual path that the particle goes through depends on the angle (the detector has cylindrical, 
but not spherical, symmetry): for each angle we have a different travel length. 
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B. dE/dX 



When a charged particle goes through the detector it deposits energy by ionization. The 
amount of energy deposited, dE/dX, is a function of fij of the particle fl9fl . We assumed 
that the charged particle goes through argonne and dE/dX has 5% resolution, which is a 
realistic value for a TPC tracking chamber. To remove muons, we use the following cut, 

fi dE/dX-dE/dX(muon) q 
D - a(dE/dX) > 6 - 

The resulting efficiency is shown as the solid line in Fig. |^a. We note a blind spot for masses 
around 150 GeV. We also note that for lower value of fij (for masses below 150 GeV) dE/dX 
will give a unique value for f3j, thus knowing the momentum we can get the stau mass. As 
the figure shows, we can apply cut 6 to cover the high mass range up to the beam energy 
where the cuts 2 and 3 are not met. 

In Fig. [|b we present, for each strategy, the minimum cross section (beforee cuts) that 
will be visible at a */s = 500 GeV linear collider with L = 50 fb _1 . Our criteria is based 
on a 3 sigma significance; namely, S > 3 with S = eery 1 L/bg, where a is the signal cross 
section (before cuts), e is the efficiency to pass the cuts and bg is the expected background 
cross section after cuts. We require a minimum of 5 signal events after cuts. 

A Cerenkov device can be used to measure [3j. With a device similar to the BABAR 
detector |2(J it is possible to reject particles with fiy > 8. Cerenkov devices, however, in 
general have a large impact on the detector design due to its requirements on space and 
photon detection. If these requirements are met, a Cerenkov device can replace TOF or 
dE/dX discussed above. 

A comment on the nature of our results is in order. We have presented a strategy 
based only on the pair production mechanism where particle identification had a relatively 
minor role to play. Nevertheless, in the models under consideration it is likely that others 
super symmetric particles will be produced and end up in stau; in such decay chains the use 
of time of flight, dE/dX and Cerenkov devices would play a critical role in identifying staus 
0. We also note that the results presented so far in Fig. |3] are rather independent of model; it 
depends only in the pair production mechanism that will exist for any non-colored charged 
long-lived particle, provided that the cross section is not too small. The efficiency to pass 
the cuts depends only on the mass of the particle 0, so that Fig. ^| can be read in a model 
independent way: we see the minimal cross section that would lead to an observable signal 
as a function of the mass of a long-lived charged particle, regardless of the nature of this 
particle. 



2 The kinematical distribution, as proposed here, would be of no use in a process other than pair 
production. 

3 There is a model dependent part in the angular distribution of the cross section that will have 
an effect in cut 1. But we note that the effect of this cut is small. One can also read Fig. ^ as the 
minimum cross section in the central region that will be observed at the NLC. 
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III. STAU DECAY LENGTH 



Up to this point we supposed that the stau does not decay within the detector; however, 
as mentioned earlier, the lifetime can be viewed as a free parameter of the model. We will 
now discuss a simple measurement of the stau lifetime using the mode e + e~ — > f + f~. The 
case of a short-lived stau is to be studied elsewhere. 

In order to ensure good measurement of momentum and dE/dX we require that each 
event should have one central track longer than 1 meter, using the mass cut 4 to select tracks 
consistent with a heavy particle. We believe that such events will be essentially background 
free: a well reconstructed track, consistent kinematically with a heavy particle, regardless 
of the decay pattern. 

One way to get the lifetime is to consider the number of events that decay before and 
after a certain length. With this method the error in the lifetime is given by: 



Or 



' VR (6) 



cr y/N]og(l + R) 

where Ni is the number of particles that decays between distance l\ and l 2 , A inf is the 
number of particles that decays after distance 1% and N = N\ + A^ n f is the total number of 
particle that decay after distance l\. In practice, l\ = lm, and 1% is chosen such that the 
statistical power is maximized, which is found to occur at R — 3.9, but not exceeding the 
outer radius of the tracking volume. When the lifetime is short enough the optimum value 
of li occurs within the tracking volume, and the error [Eq. ([])] is given by On the other 
hand, a large lifetime may lead to the optimum Z 2 outside of the tracking volume in which 
case the error is approximately In both cases the relevant quantity to get a precise 

measurement of the lifetime is the number of particles that decay inside the detector. 

In Fig. [|, the solid lines represent the contour plot of the lifetime resolution as a function 
of the stau mass and lifetime where minimum track length of 1 m is required. In many 
models, stau can be identified without using its own track. In order to take these cases 
into account, we also show the lifetime resolution for l± = 1 cm (the dashed lines). We see 
from this plot that in the optimistic case we would be able to measure the lifetime with a 
10% precision as long as cr is lower (bigger) than ~ 40 m (~ 0.01 m), for masses up to 
200 GeV. A comparison with the CERN LHC is in order O]. In this study Ambrosanio et 
al. show the capability of the CERN LHC to measure the stau lifetime in several models 
within the GMSB context. Their numbers are parameter space dependent as all possible 
reactions (that ended in staus at the end of the decay chain) are used to extract the stau 
lifetime. It is of fundamental importance for the LHC to consider all possible reaction in 
order to get enough statistics. Their analysis shows that for models that gives M stau ~ 100 
GeV they are able to get 10% precision (in a somewhat optimistic case) for 1 < cr(m) < 50. 
Our analysis involves only the stau pair production, being essentially model independent Q. 



4 We note that we used an almost pure right-handed stau that indeed gives a somewhat lower 
cross section than a left-handed stau. 
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For M stau ~ 100 GeV we are able to get 0.0013 < cr(m) < 50 (0.13 < cr(m) < 25) for the 
minimal track length of 1 cm (1 m) requirement. We also note that we are using a somewhat 
conservative integral luminosity of 50 fb _1 . 

Also in Fig. |] is shown in dotted lines a 90% confidence level upper limit for cr, assuming 
that stau decays outside the detector. From Eq. (Q), we can see that this value corresponds 
to a lower bound yF ~lx 10 7 GeV, which is slightly above cosmologically preferred value. 

IV. CONCLUSIONS 

In many supersymmetry scenarios a charged long-lived particle is predicted. For instance, 
in the GMSB scenario such particle would be a stau in a large portion of parameter space. We 
have studied the long-lived stau pair production in a linear collider at ■y/s = 500 GeV. The 
linear collider will be able to study the stau pair where no other particles are produced in the 
event (except for beamsstrahlung, etc.), as opposed to the LHC where all super symmetric 
reactions should be taken into account [0,0 , being in principle a good place to extract the 
parameters of the model. 

We presented in a model independent way the minimal cross section for the pair produc- 
tion of stable non-strongly interacting charged particle that will be observed in the NLC, 
considering just momentum measurement as well as particle identification devices. When 
the predicted cross section for the stau pair production in the GMSB model is considered 
it is shown that the NLC will be able to detect such reactions for stau masses up to 85% 
of the beam energy, with just momentum measurement. Particle identification devices will 
extend the mass range to essentially the full beam energy. Moreover, particle identification 
devices will provide a sample of events essentially background free. 

We also presented a way to extract the stau lifetime using the predicted cross section 
and a luminosity of 50 fb _1 . The method presented is very straightforward and does not 
depend in any other parameter of the model. For the range of masses that is possible to 
probe in a 500 GeV linear collider, the precision obtained is, in general, better then it will 
be possible to have from the LHC and it is mostly model independent. 
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FIG. 1. Stau pair production cross section normalized by muon pair as a function of stau mass 
in y/s = 500 GeV linear collider. The full (dotted) line is for the right- (left-) handed stau pair. 
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FIG. 2. The mass square distribution for the muon and for several values of stau masses. 
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FIG. 3. a) The efficiency for the signal. The dotted line is the efficiency for just kinematics 
cuts (cuts 1-4), the dashed line is for the time of flight cut (cuts 1 and 5) and the solid line is for 
dE/dX cuts (cuts 1 and 6). b) The reach in cross section for a linear collider with luminosity of 50 
fb _1 . The dotted line is using just kinematics, the dashed line is using time of flight and the solid 
line is using dE/dX. 
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FIG. 4. Contours of constant error a CT on the measurement of or. The solid lines stand for 
the case l± = 1 m and dashed lines are for l\ = 0.01 m. The dotted line on the top of the figure 
indicates at 90% CL the case stau is stable up to the detector. 
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